arXiv: 1502.01819v2 [physics.optics] 23 Mar 2015 


Coupled-mode induced transparency in aerostatically-tuned 
microbubble whispering gallery resonators 

Yong Yang,^'0 Sunny Saurabh,^ Jonathan Ward,^ and Sfle Nic Chormaic^ 

^Light-Matter Interactions Unit, Okinawa Institute of Science and Technology, Onna, Okinawa 904-0495, Japan 


compiled: March 24, 2015 


Coupled-mode induced transparency is realized in a single microbubble whispering gallery mode resonator. 
Using aerostatic tuning, we find that the pressure induced shifting rates are different for different radial order 
modes. A finite element simulation considering both the strain and stress effects shows a GHz/bar difference 
and this is confirmed by experiments. A transparency spectrum is obtained when a first order mode shifts 
across a higher order mode through precise pressure tuning. The resulting lineshapes are fitted with the theory. 
This work lays a foundation for future applications in microbubble sensing. 
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In the last few decades, a category of optical micro¬ 
cavities, termed whispering gallery mode (WGM) res¬ 
onators (WGRs), has been developed and has gener¬ 
ated interest in a wide range of fields. These resonators 
are characterized by having ultrahigh quality (Q) factors 
and small optical mode volumes, thereby making them 
advantageous for studies in areas such as cavity QED 
0,i, nonlinear optics Q, optomechanics and high 
sensitivity sensing Within the WGR family, the mi¬ 
crobubble resonator [g is a relatively recently-developed 
type. Experimentally, microbubbles can be made with 
single 0 or double inputs @ so that different materials 
can be injected into the hollow core. A subtle change 
in the material can cause a large frequency shift to the 
resonator’s modes, which can be used for sensin g. Exten- 
sive work has already be done in the field of bio and 
chemical sensing Q using WGMs in hollow optical res¬ 
onators. Various microbubble sensing applications have 
been studied, such as refractive index [l0|, temperature 
[HI, and optomechanical sensing [l2|. Aside from sens¬ 
ing, lasing in hollow cavities has been achieved [T^. H^. 

Another potential application of the microbubble is 
that of coupled-mode induced transparency. Electro¬ 
magnetic induced transparency (EIT)-like and Fano-like 
lineshapes in the transmission spectrum of a WGM res¬ 
onator coupled to a waveguide, such as a tapered op¬ 
tical fiber, have been observed in recent years. These 
phenomena are very important for applications such as 
slow light [T^ and optical switching [I^. To date, EIT- 
and Fano-like lineshapes have been seen in a single mi¬ 
crotoroid [13, [13 , a single microsphere (l9j | , and in di¬ 
rectly [23 or indirectly coupled [2l[ resonators. For a 
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single resonator, a low and a high Q mode should be 
found within a very small frequency separation. When 
the two modes are tuned into resonance, an EIT-like 
lineshape occurs. 

In this paper we will show that resonance can be 
achieved by pressure tuning a single microbubble res¬ 
onator [23, so that two modes with different mode or¬ 
ders experience different shift rates. Aerostatic tuning 
of the modes has some advantages over other methods. 
For example, it is precise and fast, with a larger tun- 
ing range than dispersion tuning in a single microsphere 
[l9|- Also, it is more convenient and repeatable than 
temperature tuning [TTI |. 

WGMs with different orders have different optical 
paths, or different effective mode radii, in the cavity. 
Therefore, modes of different orders show different FSRs 
and responses to external perturbations, such as aero¬ 
static pressure. These conditions are necessary for the 
EIT phenomena reported in this paper, because they 
lead to modes with close frequencies and varying aero¬ 
static tuning rates. Gonsider a microbubble of inner and 
outer radius, b and a, respectively. A microbubble with 
a thick wall, i.e., (a — b) >> A, can be treated equivalent 
to a microsphere, at least for the radial modes <7 = 1 and 
2 when the resonances are determined as an asymptotic 
formula [23 ■ The sensitivity is given as 
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where n is the refractive index of the microbubble, 
= Z -|- 1/2 and I is the angular momentum of the reso¬ 
nance. Aiq is the root of the Airy function, Ai{—z). 
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Note that in high order terms in the asymptotic ex¬ 
pression were omitted so that dX/X ~ da/a + dn/n. For 
the two modes, 9 = 1 and 2, with the same Aq, the shift 
rates differ by 


^ ^dX^ _ Xq dn^ 
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Here, Aii < Ai 2 and, for a fixed Aq, > h in general, 
so that Ail < A^2 For a silica microbubble, 

n = 1.4457. Eq. [5] is always positive, therefore low 
order modes shift faster than higher ones. The larger 
the value of dn/n, the larger the shift rates between the 
modes. From [22j| . dn/n is proportional to a geometrical 
parameter, x = /{a^ — b^), and we can write 

dn _ 3{p^b^ - poa^) ^ _ 3xC'(p^ - Po) 

/ Q 7 Q \ ^ ~ 7 

n n[a'^ — b'^) n 

where C = 4 x 10“^^ m^/N is the elasto-optic constant 
of silica, and pi and Po are the internal and external pres¬ 
sures. Therefore, for a thinner wall, the effect should be 
more obvious. However, the earlier assumption that the 
mode distribution is equivalent to that of a microsphere 
is not valid for thin-walled microbubbles; in the follow¬ 
ing discussion, we adapt a finite element method [ 1 ^ 
to precisely simulate the microbubble modes. The fact 
that the strain contribution will lead to a much larger 
shift rate difference in thin-walled bubbles should also 
be considered and we will discuss this later. 

A microbubble was made using the method described 
in m with one input sealed using a CO 2 laser. The ra¬ 
dius and wall thickness were measured as a = 45 ^m and 
t = a — b = 1.26 pm. The WGM shifts were recorded 
using a taper-coupled setup illustrated in Fig. [IJa) and 
an electronic sensor was calibrated for monitoring the 
internal pressure. The taper was made by heating-and- 
pulling standard 1550 nm single-mode fiber down to a 
waist of less than 1.2 pm to ensure single mode oper¬ 
ation. To eliminate the influence of vibrations, the mi¬ 
crobubble was in contact coupling with the tapered fiber 
and the whole system was placed in an enclosure. Two 
modes with different Q factors were chosen for tracking. 
It is assumed that the low Q mode was q = 2 and the 
high Q mode was 9 = 1. Resonant frequencies at differ¬ 
ent internal pressures are plotted in Fig. [2ja). The linear 
fits yield shift rates of 4.27 GHz/bar and 2.60 GHz/bar 
for the two modes, resulting in a positive rate difference 
of 1.67 GHz/bar. 

To simulate the experimental results, the frequency 
shifts due to stress and strain were calculated separately. 
By setting Pi — po in Eq. |3]to 0, 1 x 10^ Pa (i.e., 10 bar) 
and 2x 10"^ Pa (i.e., 20 bar) we numerically calculated the 
frequency shifts for the 9=1 and 2 modes (Fig[TDb)) in 
a bubble with a diameter of 90 pm and a wall thickness 
varying from 1 pm to 1.6 pm. The stress effect, dn/dp^ 
causes a shift rate difference in the MHz/bar scale only. 



Fig. 1. a) Schematic of the experiment. TLS: tunable laser 
sonrce; DSO: digital oscilloscope; PD: photon detector. P is 
an electrical pressure sensor. The inset is the photo of the 
microbubble used in the experiments, with a diameter of 90 
^m and wall thickness of 1.26 pm. b) The simulated EM held 
distribution for the q = 1 and 2 modes (from left to right) 
using a hnite element method. 


Strain affects the WGM resonance by changing the 
size of the microbubble as da/dp. FEM simulations were 
used to calculate the shifts for the 9 = 1 and 2 modes 
for a given size expansion (i.e., 0.01 pm increase in a). 
Next, the relative change in a, given by 
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was used to determine the required pressure for the 
same expansion. This was repeated for different wall 
thicknesses. The shear and bulk moduli in EqH] are 
G = 31 X 10® Pa and AT = 41 x 10® Pa for the silica 
capillary used. As a summary. Fig [5Kb) illustrates the 
shift rate differences between the two modes originat¬ 
ing from the strain and the stress, thus giving the total 
rate difference. For thin-walled microbubbles the strain- 
induced rate difference is in the GHz/bar range. The 
shift rate difference for the 9 = 1 and 2 modes for a wall 
thickness of 1.3 pm is calculated to be 1.90 GHz/bar, 
which is quite close to the measured value from Fig. 
Ha). 

In experiments, the linewidth of the low Q mode is 
about 2 GHz. We could always find a high Q mode 
(linewidth of MHz) in the vicinity of the low Q one. 
With only about 1 bar of input pressure, the high Q 
mode (A) can be tuned across the low Q mode (B) 
transmission dip. Fig. shows an example of such a 
crossing, where an EIT-like shape appears as a result of 
tunin g th e two modes into resonance. 

In [l8|, the appearance of EIT and Fano shapes in 
a single microresonator was explained by coupled-mode 
theory as follows: 


dUA 

—— = -ILOAttA - 
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A and B are considered to be two quasi TE (TM) WGMs 
indirectly coupled through the taper mode a™. xa,b 
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Fig. 2. (a) Experimentally measured pressure tuning rate 
of a high Q mode (squares) and a low Q mode (spades) in 
a microbubble. The solid red and blue lines are linear fits 
to the data. (b)Frequency shift rate differences of g = 1,2 
modes for 1-1.6 /rm thicknesses. The black squares are due 
to the strain effect and red circles are total rate differences. 
Lines are only for guiding the eye. 


are the intrinsic rates and 0 is the polarization ori¬ 
entation of the taper mode a™, see Fig. [3Ka). The 
coupling strength is determined by the external cou¬ 
pling rates, ka and kb, and their polarization angles, 
if, through g = {^Jhaub cosip)/2. Intracavity EM 
fields of the A, B modes, i.e. a a and ob, are cou¬ 
pled back to the taper fiber mode, with a different 
input/output relationship for two orthogonal polariza¬ 
tions as a™* = a*" cos 9 -I- ^/nAaA + \/kb cos ipa b and 
a°"‘ = a™ sin 9 + sin ipa b ■ The transmission spec¬ 

trum is an intensity superposition of these to orthogonal 
polarizations: 


T = 


( 6 ) 


From Eq. [5l aA.B can be deduced when in steady state, 
i.e. da/dt = 0. By substituting into Eq. [51 the transmis¬ 
sion can be calculated for known [ka, kb,^a,^b, P>,9] 
and detuning Aui = lub — to a- Thus EIT and Fano-like 
lineshapes can be manipulated by either changing the 
detuning or polarization [T^. 11^. 

In contrast to earlier work [19l| , where the detuning of 
the two modes is modified by varying the conpling gap 
and, subsequently, by dispersion, our system remains 
in contact with the coupling fiber. The pressure tun¬ 
ing only changes the detuning of the modes with a mnch 
larger range, while keeping other parameters unchanged. 
Since the rate difference has been determined in the pre¬ 
vious discussion, the detuning can be precisely known if 
the pressure is also known. In practice, the experimental 



Fig. 3. (a) A schematic illustration of coupled mode model of 
a microbubble coupled with taper, (b) Fitted EIT-like shape 
with pressure tuning. The red curves are theoretical fits with 
fitting parameters [ka, 7a, ks, 7s] = [0.1,0.08,1.0,2.2] GHz, 
ip = 7r/10 and 6 = 7r/6. From top to bottom the detunings 
are 0.25, 0.1, 0, -0.1 and -0.25 GHz. 


data was fitted with the detuning ranging from 250 MHz 
to -250 MHz and using one unchanged set of parameters, 
9], which are given in Fig. [31 Note there 
is still a small discrepancy when adding more pressure 
(i.e. detuning is more negative). A possible explana¬ 
tion is that silica is a birefringent material, so pressuriz¬ 
ing the microbubble may rotate the polarization of the 
WGMs, thereby distorting the lineshape from EIT-like 
to Fano-like. 

In summary, we have shown that a controllable EIT- 
like shape can be easily realized in a single aerostatically- 
tuned microbubble. We proved, both theoretically and 
experimentally, that, for different radial modes, the pres¬ 
sure tuning rates are different due to the size and refrac¬ 
tive index change induced by strain and stress. The ex¬ 
periments show that aerostatic tuning of the microbub¬ 
ble might be a simplified and precise way to implement 
controllable EIT and Fano lineshapes. Furthermore, 
by having a Fano-like lineshape due to the coupling of 
two modes, the sensitivity for pressure sensing in a mi¬ 
crobubble can be enhanced, as predicted in |I8| . 
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